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Abstract 



H I The role of R in the solutions of the Ginsparg- Wilson relation is 

discussed. 
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In this note, we would like to clarify some seemingly subtle issues pertaining 
to the role of R in the solutions of the Ginsparg- Wilson relation [|l| 

^75 + I5D = 2D^^RD . (1) 

Here i? is a positive definite hermitian operator which is local in the position 
space and trivial in the Dirac space. Since we can sandwich (|T|) by a left 
multiplier yR and a right multiplier yR on both sides of (|I|), and define 
D' = VRDVR, then D' satisfies 

D'75 + 75l^' = 2D'75^' (2) 

which is in the same form of Eq. (|l]) with R = 1. This seems to suggest 
that one can set i? = I in the GW relation (0) and completely ignore the R 
dependence in the general solution of the GW relation. However, as we will 
see, if R is set to be the identity operator from the beginning, then some of 
the salient features of the general solution may be easily overlooked. 

First, let us review some basics of the GW relation. The general solution 
to the GW relation (|^) can be written as [^ Q 

D = D,{1 + RD,)-' = (I + D.Ry'D, (3) 

where D^ is any chirally symmetric Dirac operator, i.e., 

D,75 + 75^c = . (4) 

In order to have D reproduce the continuum physics, Dc is required to satisfy 
the necessary physical constraints [Q. The general solution of D^ has been 
investigated in ref. @]. Conversely, for any D satisfying the GW relation (0), 
there exists the chirally symmetric Dc 

Dc = D{1- RD)-^ = {1-DR)-^D . (5) 

We usually require that D also satisfies the hermiticity condition^ 

D^ = 75D75 • (6) 

Then Dc also satisfies the hermiticity condition Dl = •y^Dc'j^, since R is her- 
mitian and commutes with 75. The hermiticity condition together with the 
chiral symmetry of Dc implies that Dc is anti-hermitian. Thus there exists one 
to one correspondence between Dc and a unitary opertor V such that 

Dc = il + V)il~Vy\ V = {Dc-l)iDc+l)-\ (7) 



^ This implies that det(£') is real and non-negative. 



where V also satisfies the hermiticity condition V"^ = 751^75. Then the general 
solution (H) can be written as H, 0] 

D = {l + V)[{l-V) + R{l + V)]-' (8) 

= [{1-V) + {1 + V)R]-\1 + V). (9) 

On the other hand, if one starts from Eq. @), then its solution is 

D' = ^il + V) (10) 

which agrees with (H) with R = 1. Then using the relation D' = yRDyR^ 
one obtains 



D 



^ (I + V^) ^ 



'2R 



2R 



;ii) 



However, ([Tl| ) is in contradiction with @ since the R dependence in ( [TlD 
can be factored out completely while that of (^ can not. So, one of them 
can not be true in general. If we take the limit R —>■ 0, then (H) gives that 
D^ {1 + V){1-V)-^ = De, but (0) implies that D ^ 00. Since D^ is well 
defined ( without poles ) in the trivial gauge sector, it follows that ( pTj ) can 
not be true in general. The fallacy in ([TTl) is due to the assumption that D' is 
independent of R. From (Rl), one can derive the following formula 
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' 'i + r- ' 



2R 



2R 



where 
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'I + V)VR+il-V)^ (l + V)VR-il-V)^ 



R{l + V) 



(l-V) 



R{l + V) 



(l-V) 



(12) 

(13) 
(14) 



which is unitary and depends on R. This shows that D' actually depends on R 
and equals to ^{l + V) rather than (|T0|). Therefore it is erroneous to write the 
general solution of the GW relation in the form of (0), in which V is replaced 
by V. Consequently, (^) may mislead one to infer that R does not play 
any significant roles in the locality oi D, in particular when R is proportional 
to the identity operator. However, R indeed plays a very important role in 
determining the locality of D. Let us consider R = rl with r > 0. When 
r ^ 0, D ^ Dc which must be nonlocal if Dc is free of species doubling and 
has the correct behavior in the classical continuum limit 0. As the value of 
r moves away from zero and goes towards a finite value, D may change from 
a non-local operator to a local operator. This has been demonstrated in ref. 



[|, 0]. Therefore, the general solution (^) of the GW relation can be regarded as 
a topologically invariant transformation ( i.e., mdex{D) = mdex{Dc) ) which 
can transform a nonlocal D^. into a local D. Conversely, the transformation 
(^ can transform a local D into the non-local D^. For a given D^, the set of 
transformations, { T{R) : D = Dc{l + RDc)~^ }, form an abelian group with 
parameter space { i? } [0. In general, for any lattice Dirac operator D ( not 
necessarily satisfying the GW relation ), we can use the topologically invariant 
transformation D' = D{1 + RD)~^ to manipulate its locality. 

The next question is whether we can gain anything ( e.g., improving the 
locality of D ) by using another functional form of R rather than the simplest 
choice R = rl. We investigate this question by numerical experiments. For 
simplicity, we consider the Neuberger-Dirac operator [Q 

Dh = l + V, V = D^{DlD^)-^'^ (15) 

where D^) is the Wilson-Dirac fermion operator with negative mass —1 

D^ = -l + i[7;.(V* + V^) - V* V^] (16) 

where V^ and V* are the forward and backward difference operators defined 
in the following, 

V/,V(x) = U^{x)i){x + jx) -i){x) 
V*f,i){x) = ij{x) - Ul{x - fi)i){x - ft) . 

The Neuberger-Dirac operator D^ satisfies the GW relation (|l|) with R = 1/2. 
Then Eq. (^) gives 

Substituting ( |T7|) into the general solution (|^), we obtain 

D = 2{1 + V)[{1-V) + 2R{1 + V)]-^ . (18) 

For a fixed gauge background, we investigate the locality of D{x, y) versus the 
functional form of R{x, y). For simplicity, we consider D in a two dimensional 
f/(l) background gauge field with non-zero topological charge, and we use the 
same notations for the background gauge field as Eqs. (7)-(ll) in ref. 0. 
Although it is impossible for us to go through all different functional forms of 
R{x,y), we can use the exponential function 

R{x,y) = r exp(— m|x — yl) (19) 

as a prototype to approximate other forms by varying the parameters r and m. 
In the limit m ^ 0, R{x,y) is nonlocal, while in the limit m -^ oo, R ^ rl 



which is the most ultralocal. Hence, by varying the value of m from to 
a ^ 1, we can cover a wide range of R{x, y) of very different behaviors. 

One of the physical quantities which are sensitive to the locality of D is 
the anomaly function 

A(x) = tr Yi^{RD){x, x) + i^{DR){x, x)] (20) 

which can serve as an indicator of the localness of D. Since the Neuberger- 
Dirac operator is topologically proper for smooth gauge backgrounds, the index 
of D in ( p!8D is equal to the background topological charge Q, 



index(D) = n_ - n+ = Q . (21) 

This implies that the sum of the anomaly function over all sites on a finite 
lattice must be equal to two times of the topological charge [0 



^A(a;) = 2(n_-n+) = 2g 






(22) 

This is true for any R since the index of D is invariant under the transformation 
(H), i.e., index(D) = index(Dc)- First we consider the gauge configuration with 
constant field tensors. If D is local, then we can deduce that An{x) is constant 
for all X. From (E2D, it follows that 



An{x) = p{x) = { (23) 

where p{x) is the Chern-Pontryagin density in continuum. Note that Eq. (^31) 
also implies that An{x) is independent of i? if Z) is local. Next we introduce 
local fluctuations to the constant background gauge field, with the topological 
charge fixed. Then we expect that (p3| ) remains valid provided that the local- 
ity of D is not destroyed by the roughness of the gauge field. Therefore, in 
general, by comparing the anomaly function An{x) at each site with the Chern- 
Pontryagin density p(x), in a prescribed background gauge field, we can reveal 
whether D is local or not in this gauge background. This provides another 
scheme to examine the locality of D rather than checking how well \D{x,y)\ 
can be fitted by an exponentially decay function. We will use both methods 
in our investigations. We define the deviation of the anomaly function as 

where the summation runs over all sites, and A^, is the total number of sites. 



r 


m 


Sd 


index (D) 


0.5 


0.5 


1.366 




0.5 


1.0 


0.5360 




0.5 


2.0 


0.2086 




0.5 


5.0 


7.209 X 10-^ 




0.1 


5.0 


6.554 X 10"^ 




0.5 


>1 


3.848 X 10--* 




1.0 


>1 


1.618 X 10-"* 





Table 1: The deviation of the chiral anomaly function 6d [ Eq- (P^) ] versus 
R{x,y) with parameters r and m defined in Eq. ([T9|) , on a 12 x 12 lattice, in 
a constant background gauge field with topological charge Q = 1. 



In Table |l|, we list the deviation of the chiral anomaly function, 6d, versus 
R{x,y) with parameters r and m defined in Eq. ([T9|) , on a 12 x 12 lattice, in 
a constant background gauge field with topological charge Q = 1. The last 
two rows with ?7i ^ 1 corresponds to R = rl, and they have the smallest 
deviations. They are both local, which can be checked explicitly by plotting 
\D{x, y)\ versus |x — 1/|, as shown in Fig. 1. For r = 1/2 ( the row on top of the 
last row ), it corresponds to the Neuberger-Dirac operator. In the first row, 
R{x,y) is nonlocal since m = 1/2 < 1. It produces a nonlocal D, as shown 
in Fig. 2, so the resulting chiral anomaly is very different from the Chern- 
Pontryagin density and thus 6d is very large. Now we increase m to 1, 2 and 
5 successively, then R and D both become more and more local, as shown in 
Fig. 3, thus 6d becomes smaller and smaller, as shown in the second, third 
and fourth rows of Table 1. These results indicate that a nonlocal R does not 
produce a local D, and a local R does not make D more local than that using 
R = rl. For m = 5.0, if we decrease r to 0.1, then D becomes very nonlocal, 
as shown in Fig. 4. Consequently, its Sd ( in the fifth row ) is about 10 times 
larger than that of r = 0.5 ( in the fourth row ). This suggests that on a finite 
lattice, r cannot be too small, otherwise D will become nonlocal. 

Our numerical results listed in Table 1 as well as those plotted in Figs. 1-4 
strongly suggest that we do not gain anything by using other functional forms 
of R{x,y) than the simplest choice R{x,y) = r 5x,y However, the value of r 
plays the important role in determining the localness of D. We have also tested 
other functional forms of R{x, y) as well as many different gauge configurations. 
The results from all these studies are consistent with the conclusion that the 
optimal choice for R is R{x, y) = r 5x,y 

Now we come to the question concerning the range of proper values of r. 
We have already known that r cannot be zero or very small, otherwise D is 
nonlocal. On the other hand, r cannot be too large, otherwise D is highly 
peaked in the diagonal elements ( i.e., Dap{x,y) ~ Daa{x,x)6ai3Sx,y ), which 



r 


Sd 
( 12x 12 ) 


( 16 X 16 ) 


Sd 
( 20 X 20 ) 


index (D) 


0.1 


6.232 X 10-^ 


1.792 X 10-^ 


4.854 X 10-^ 




0.2 


4.818 X 10-^ 


6.229 X 10-* 


7.750 X lO-'' 




0.5 


3.848 X 10-"* 


2.434 X 10"^ 


1.559 X 10-" 




0.8 


1.698 X 10-" 


8.268 X 10-" 


4.976 X 10-^ 




1.0 


1.618 X 10-"* 


6.596 X 10-" 


3.432 X 10-^ 




1.2 


3.448 X 10"'' 


1.318 X 10-^ 


5.308 X 10-' 




1.5 


1.688 X 10-^ 


1.243 X 10-" 


8.783 X 10-" 




2.0 


9.940 X 10-^ 


1.569 X 10-^ 


2.348 X 10-" 




5.0 


0.1599 


9.314 X 10-^ 


5.088 X 10-" 





The deviation of the chiral anomaly function, (5b [ Eq. 



versus 



Table 2: 

The background gauge field has constant field tensors with topological charge 
Q = 1. The index of D is always equal to one in each case. 



6x,y, for lattice sizes, 12 x 12, 16 x 16 and 20 x 20 respectively. 



is unphysical since it does not respond properly to the background gauge field 
( e.g., the chiral anomaly is incorrect even though the index of D is equal to 
the background topological charge ). In Table 0, we list the deviation of the 
chiral anomaly function, 6d, versus R{x,y) = r 6x,y, on a 12 x 12 lattice ( the 
second coluum ), in a constant background gauge field with topological charge 
Q = 1. We see that the proper values of r are approximately in the range 
0.5 ~ 1.2, where D can reproduce the continuum chiral anomaly precisely. 
Next we investigate how the lattice size affects the range of proper values of r. 
The results of Sn for lattice sizes 16 x 16 and 20 x 20 are listed in the third and 
the fourth coluums in Table ^ They clearly show that the lower bound of r 
can be pushed to a smaller value, ~ 0.2, when the size of the lattice is increased 
to 20 X 20. Therefore, it suggests that the chiral limit ( r — > and D ^ Dc) 
can be approached by decreasing the value of r while increasing the size of the 
lattice, at finite lattice spacing. This provides a nonperturbative definition of 
the chiral limit for any D of the general solution (^) with D^ satisfying the 
necessary physical requirements 0. 

It is evident that the range of proper values of r also depends on the back- 
ground gauge configuration. However, we suspect that when the background 
gauge configuration becomes very rough, there may not exist any values of r 
such that the chiral anomaly function is in good agreement with the Chern- 
Pontryagin density. We intend to return to this question in a later publication. 

In summary, we have clarified the role of R in the general solution (^ of the 
Ginsparg- Wilson relation. It provides a topologically invariant transformation 
which transforms the chirally symmetric and nonlocal Dc into a local D which 
satisfies the GW relation, the exact chiral symmetry on the lattice. Having R 
local in the position space is a necessary condition to ensure the absence of ad- 



ditive mass renormalization in the fermion propagator, as well as to produce a 
local -D, which is vital for obtaining the correct chiral anomaly. Our numerical 
results strongly suggest that the optimal form of R is R{x,y) = r d^^y The 
range of proper values of r depends on the background gauge configuration 
as well as the size of the lattice, L = Na. In the limit A^ — > oo, for smooth 
gauge backgrounds, the lower bound of proper values of r goes to zero, thus 
the chiral limit ( r ^ and D ^ Dc) can be approached nonperturbatively 
at finite lattice spacing. 
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Figure 1: One of the Dirac components of D{x,0), \Dii{x,0)\, is plotted 
as a function of |a;| for R{x,y) = Sxy The lattice is 12 x 12 with periodic 
boundary conditions. The constant background gauge field has topological 
charge Q = 1. All data points at the same distance |x| from the origin have 
been averaged. The solid line is an exponential fit to the data points. The same 
decay constant also fits very well for all other Dirac components of D{x, y) and 
for any reference point y. 
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Figure 2: One of the Dirac components of D{x, 0) is plotted as a function of 
|a;| with R{x,y) = 0.5exp(— 0.5|a; — y\). Other descriptions are the same as 
Fig. 1. The non-localness of D is shown clearly. 
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Dii(x,0) 



R(x,y) = 0.5 exp( -5.0 | x - y | ) 



ly = 2.034exp{-2.143x; 
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Figure 3: One of the Dirac components of D{x, 0) is plotted as a function of 
|a;| with R{x,y) = 0.5exp(— 5.0|a; — y\). Other descriptions are the same as 
Fig. 1. 
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D„(x,0) 



R(x,y) = 0.1 exp( -5.0 | x - y | ) 
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Figure 4: One of the Dirac components of D{x, 0) is plotted as a function of 
|a;| with R{x,y) = 0.1 exp(— 5.0|a; — y\). Other descriptions are the same as 
Fig. 1. The non-localness of D is shown clearly. 
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